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We report the measurement of charged D* mesons in inclusive jets produced in proton-proton 
collisions at a center of mass energy yfs = 200 GeV with the STAR experiment at RHIC. For D* 
mesons with fractional momenta 0.2 < z < 0.5 in inclusive jets with 11.5 GeV mean transverse 
energy, the production rate is found to be N(D* + + D* _ )/W(jet) = 0.015± 0.008(stat) ± 0.007(sys). 
This rate is consistent with perturbative QCD evaluation of gluon splitting into a pair of charm 
quarks and subsequent hadronization. 

PACS numbers: 13.85.Ni, 13.87.Fh, 13.25.Ft 



Charm and bottom quarks can probe the par- 
tonic matter produced in heavy-ion collisions [ij 
and the nucleon spin structure in polarized proton- 
proton collisions [2J. Their production mechanism 
is, therefore, of considerable interest at the Rela- 
tivistic Heavy Ion Collider (RHIC). Studies of the 
Z?* ± -meson content injets by the ALEPH, L3 and 
OPAL Collaborations [3j show that the production 
from Z° decays in e + + e~ collisions is dominated by 
D* mesons that carry large fractions of the jet mo- 
menta, consistent with the jets being produced from 
primary c (anti-) quarks. The E531 and NOMAD 
Collaborations observed events with large D* + mo- 
mentum fractions in neutrino charged-current inter- 
actions 4] . In p + p collisions at yfs — 630 GeV and 
1.8 TeV, the UAl and CDF Collaborations have ob- 
served D*^ mesons in jets with transverse energies 
larger than 40 GeV H Q . Their fractional momenta 
are found smaller, consistent with a different produc- 
tion mechanism in which the D* mesons originate 
from gluon splitting into cc pairs (g — > cc in the ini- 
tial or final parton shower, with neither of the quarks 
from the ccpair participating in the hard QCD in- 
teraction) [7| . The multiplicity of such heavy quark 
pairs in gluon jets is calculable in perturbative QCD 
(pQCD), and the leadingnonperturbative correction 
is believed to be small Q. At yfs = 200 GeV, the 
RHIC energy, heavy quarks can still be produced 
via gluon splitting. Perturbative QCD suggests that 
these contributions are small, and that the majority 
of the heavy quarks originate from gluon-gluon fu- 
sion These expectations, however, have not 
until now been confronted with data at RHIC. 

In this paper we present the first measurement of 
charged D* mesons in inclusive jets produced in p + 
p collisions at a center of mass energy y/s = 200 GeV 
at RHIC. The data were recorded in the year 2005 
with the Solenoidal Tracker At RHIC (STAR) [ll[ 
and amount to an integrated luminosity of 2 pb _1 . 
The main subsystems used in the measurement were 
the Time Projection Chamber (TPC) and the Bar- 
rel Electro-magnetic Calorimeter (BEMC), both lo- 
cated in a 0.5 T solenoidal magnetic field. The TPC 
provided tracking for charged particles with pseudo- 
rapidities 1 77 1 < 1.3 for all azimuthal angles 4>. The 
BEMC provided triggering and was used to measure 
photons and electrons. In 2005 it covered < rj < 1 



in pseudo-rapidity and 2-7T in azimuth. Events used 
in this analysis were required to satisfy both a min- 
imum bias trigger condition and a jet patch (JP) 
trigger condition. The minimum bias trigger was 
defined as a coincidence between Beam-Beam Coun- 
ters (BBC) on either side of the interaction region, 
and the JP trigger, used also in Ref. [12j, required 
the sum of transverse energies deposited in one of 
six Arj x Acf> =1x1 patches of BEMC towers to be 
above a threshold of 6.5 GeV. 

The charged D* candidates were identified 
through the decay sequence D* + — > -D°7r+, D° — ► 
K~n + and its charge conjugate. The D* decay has 
a small Q-value. The D° thus carries most of the 
D* momentum and the pion from the D* decay, 
denoted by ir s , is soft. In the following, we will 
use D* to denote D* + + D*~ , and D° to denote 
D° + D° unless specified otherwise. The enhance- 
ment in the distribution of the invariant mass differ- 
ence AM = M(KTn ± Trf) - Af(ivT T 7r ± ) is used to 
determine the D* yield [lj|. The candidate daugh- 
ter kaons and pions were tracked with the TPC and, 
where possible, identified using the agreement of 
the measured and expected ionization energy loss 
(dE/dx) in the TPC to within two standard de- 
viations. The reconstructed tracks were required 
to have transverse momenta pr > 0.2 GeV/c and 
pseudorapidities \ij\ < 1. Only those events whose 
reconstructed primary interaction vertices were on 
the beam axis within 100 cm from the TPC center 
were retained. A mass interval 1.82 < M(K =F 7r ±) < 
1.90 GeV/c 2 was used to select D° candidates, con- 
sistent with the D° mass [l4| and the experiment 
invariant mass resolution. About 90% of D° signals 
are within this mass interval. Combinatorial back- 
ground was suppressed using the low Q- value of the 
D* decay by requiring the ratio, r, of the transverse 
momenta of the D° candidates and the soft pions 
to be 10 < r < 20. In addition, the decay angle of 
the kaon in the K^tt^ rest frame, 9* , was restricted 
by requiring cos(#*) < 0.8 to remove near-collinear 
combinatorial background from jet fragmentation. 

Figure [T] (a) shows the spectrum of the invariant 
mass difference AM = M (K^ irf 
The "right sign" combinations K^n^nf 1 were used 
in obtaining the D* candidates, while the dou- 
bly Cabbibo-suppressed "wrong sign" combinations 
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X ± 7r^7r^ = were used as a measure of combinatorial 
background. The "wrong sign" distribution in Fig. Q] 
(a) was superimposed directly on the "right sign" 
distribution, that is, without applying a relative nor- 
malization. The hatched sidebands in the K^ir - 
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FIG. 1: (a) The observed distribution of the invariant 
mass difference AM = M^K^i^irf) - M(JC =F 7r ± ) in 
p + p collisions at *Js = 200 GeV. The crosses show sig- 
nal and background, and the histograms show two eval- 
uations of the background, discussed in the text, (b) 
The invariant mass distribution of the K^w pairs for 
the events with an additional soft pion in the mass re- 
gion 143.5 < AM < 146.5 MeV/c 2 . The cross-hatched 
area depicts the D° mass interval used in selecting the 
D* candidates and the hatched areas are used in con- 
structing the background, (c) The transverse momen- 
tum distribution of the D* candidates after background 
subtraction. No corrections were applied for efficiency 
and acceptance. The lower pr bound in this spectrum 
results from kinematic selection criteria applied to the 
D* decay daughters. 



mass spectrum of Fig.[T](b) were used in an alterna- 
tive measure of the D* combinatorial background. 
In the if T 7r ± -spectrum, a subsample of all K^ir^ 1 
candidates was chosen by requiring the event to con- 
tain an additional soft pion of the right charge sign 
resulting in 143.5 < AM < 146.5 MeV/c 2 . This 
mass range contained about 95% of the D* signal. 
The sample of all K^tt^ candidates has considerably 
larger background and does not show a significant 
D° signal. The crosshatched area underneath the 
D° peak indicates the mass interval used in the re- 
construction of D* candidates. The sideband-based 
distribution of combinatorial background was nor- 
malized to the D* signal candidate distribution in 
the mass region, 150 < AM < 175 MeV/c 2 , well 
away from the D* signal. The combinatorial back- 
ground distributions from the two methods are in 
good agreement. The signal above background in 
Fig. [1] (a) has ~ 4<r significance and corresponds to 
180 ± 45 D* counts. The D*+ and D*~ yields of 
96 ± 32 and 84 ± 33 counts are equal to within their 
statistical uncertainties, as expected at this level of 
precision [l5[. The sideband-based background dis- 
tribution was used in extracting these yields since it 
results in better precision. The difference with the 
wrong-sign results was used in assessing systematic 
uncertainties of the measurement. Raw pt distribu- 
tion is shown in Fig. Q] (c). 

Jets were reconstructed using a mid-point cone 
algorithm [l6| which clusters reconstructed TPC 
tracks and BEMC energy deposited within a cone 
in r\ and <b of radius r cono = 0.4, as described in 
Refs. [TU, [TtJ • Events with reconstructed primary 
interaction vertex positions on the beam axis within 
100 cm of the TPC center were kept for further anal- 
ysis. Jets were required to have pt > 8 GeV/c, 
< r] < 1, and an electro-magnetic fraction of the 
jet transverse energy within 0.1 and 0.9 to reduce the 
effects of event pile- up and beam background (T^ |. 
A sample of 1.7 x 10 6 jets that pointed to a triggered 
jet patch was retained. 

Figure [2] shows the distribution of the D* - 
candidate azimuthal angle with respect to the re- 
constructed jet axis. Background was subtracted 
using the sideband method. The distribution was 
corrected for the D* reconstruction efficiency and ac- 
ceptance obtained from PYTHlA-based (v 6.205 [l8| 
'CDF TuneA' settings [lj|) Monte Carlo simulations 
passed through GEANT-based [20[ STAR detector re- 
sponse simulation. The same simulation setup has 
been used in Refs. [H|,[l7j and provides an adequate 
description of the inclusive jet data. The determi- 
nation of the D* reconstruction efficiency took into 
account different configurations where some or all 
of the D* decay daughters were part of the recon- 
structed jet and also different intervals for the jet 
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FIG. 2: The distribution of the D* azimuthal angle 
with respect to the reconstructed jet axis from p + p 
JP triggered data. The distribution has been corrected 
for the D* reconstruction efficiency. The curve is a two- 
Gaussian fit to the data points. 



and D* momenta were taken into account. The indi- 
cated uncertainties are the quadratic sum of the sta- 
tistical uncertainties in the data and in the Monte- 
Carlo simulation. The solid line is a two-Gaussian 
fit to the data points. A clear correlation is observed 
at the near side as expected. The away side correla- 
tion is limited by statistics. In the following we will 
focus on the near side correlation to investigate the 
production of charm in jets. 

For near side D* candidates, the fragmentation 
variable z = pu (K^^irf') / E- ]ct was computed, 
where p^K^n^nf) is the K^tt^tt^ 1 momentum 
projection on the jet axis, and £j e t is the jet to- 
tal energy. The reconstructed jet transverse energy 
is on average about 20% larger than the generated 
jet transverse energy, mostly because of the sharply 
decreasing jet yield with increasing transverse en- 
ergy and the jet transverse energy resolution. The 
resolution has been studied by Monte-Carlo simula- 
tion and by using transverse momentum balance in 
a sample of dijet events 0, The reconstructed 
jet transverse energy has been corrected and residual 
effects are accounted for in the systematic uncertain- 
ties. The distribution of z, after this correction, is 
shown in Fig. E^a) . The signal in this spectrum cor- 
responds to 72 ± 25 counts. The uncertainties rep- 
resented by the bars are statistical and the brackets 
indicate the contribution caused by combinatorial 
background subtraction. No corrections were made 
here for trigger effects and reconstruction efficiency. 
The average D* p T is ~3 GeV/c for 0.2 < z < 0.5, 
and ~6 GeV/c for z > 0.5. The average D* recon- 
struction efficiency from simulation, shown in Fig. [3] 
(b), is found to increase with increasing z. The trig- 
ger efficiency largely cancels in the measurement of 



the production ratio N(D*)/N(jet) of interest here. 
However, the JP trigger condition preferentially se- 
lects jets with large electromagnetic energy. It thus 
disfavors jets containing the hadronic decay prod- 
ucts of the D* mesons, in particular for high z. The 
effects of this trigger bias were studied by compar- 
ing the simulated jet yields with and without the JP 
trigger condition. Their ratio is found constant be- 
low z ~ 0.5 and decreases rapidly for larger z, as 
expected. The green band in Fig. [3] (a) was ob- 
tained by simulating only the direct charm flavor 
creation processes, gg — » cc and qq — > cc, in PYTHIA 
and passing the results through the STAR detector 
response simulation. The simulated data were an- 
alyzed in the same way as the real data and were 
normalized using the measured total charm produc- 
tion cross section [2lj. Only a small faction of the 
generated events containing D* mesons with z > 
0.5 satisfies the JP trigger condition. To within the 
large uncertainties good agreement is found with the 
D* data at high z, where the production of charmed 
hadrons is expected to be dominated by charm quark 
fragmentation. The excess observed in the data at 
smaller z can be ascribed to production processes 
that are not included in the simulation, such as gluon 
splitting. 

The ratio N(D*)/N (jet) was determined for the 
region 0.2 < z < 0.5. For z < 0.2, the D* re- 
construction efficiency is low, and for z > 0.5, the 
JP trigger is strongly biased against jets with D* 
mesons that decay into charged hadrons. After cor- 
recting for the D* reconstruction efficiency, shown in 
Fig. [3] (b), and the decay branching ratio of (67.7 ± 
0.5)% for D*+ -> D°tt+ and of (3.89 ± 0.05)% for 
D° -> K-tt + [13, we obtain N(D*+ + D*~)/NQet) 
= 0.015 ± 0.008(stat) ± 0.007(sys) for 0.2 < z < 0.5 
and a mean jet transverse energy of (£t) = 11.5 
GeV. The estimated statistical uncertainty includes 
the statistical uncertainty in the simulations. The 
main contributions to the systematic uncertainty are 
estimated to originate from the jet definition and se- 
lection (~35%), from trigger bias (~18%), from D* 
combinatorial background (~10%), and from the D* 
reconstruction efficiency (~10%). The uncertain- 
ties associated with the jet definition and selection 
were estimated by varying the accepted primary ver- 
tex range, the jet 77 range, and the criteria used to 
reduce the effects of event pile-up and beam back- 
ground. The effects from trigger bias were assessed 
by Monte-Carlo simulation. The size of the back- 
ground uncertainty was estimated by comparing the 
results obtained with the different background sub- 
traction methods. The uncertainty in the D* re- 
construction efficiency was estimated by varying the 
daughter particle track quality criteria. The contri- 
butions were combined in quadrature to obtain the 
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FIG. 3: (a): The distribution of the D* longitudinal 
momentum fraction z in jets from JP triggered data. 
The size of the statistical uncertainties are indicated by 
the bars and the size of the background-subtraction sys- 
tematic by the brackets. No corrections were applied for 
trigger effects and D* reconstruction efficiency; however, 
the observed jet momenta and hence z were corrected 
for the detector response. The data at large z are com- 
pared with a Monte Carlo simulation of charm creation 
through gg/qq — > eg. (b): The average D* reconstruc- 
tion efficiency versus z. 



total systematic uncertainty estimate. 

To estimate the rate of gluon splitting into charm 
pairs, Rg^ C c, from the ratio N(D*)/N(jet) one 
needs to correct for the unmeasured z region, the 
fraction of charm quarks that fragment into D* , and 
the fraction of gluon jets in the sample. The fraction 
of gluon jets in the data was estimated to be 60% 
from PYTHIA simulations and from next-to-leading 
order pQCD evaluation [22|. A 10% uncertainty 
is included as a systematic contribution in R g ^ cS . 
The c — > D* + and c — > D*~ fraction is taken to be 
(22.4 ± 2.8)% [lj]. This is smaller than the value of 
3/8 estimated in the earlier publications by UA1 Q 
and CDF @. By using the leading-order pQCD eval- 
uation of gluon splitting [8j, we estimate that the 
measured ratio N(D*)/N(jet) for 0.2 < z < 0.5 cap- 
tures (53±5)% of R g ^ c5 at (E T ) = 11.5 GeV and the 
dominant part of the remainder resides at smaller z. 
This percentage was then used to extrapolate over 
the unmeasured z region. Our result for R g ^ cS is 
shown in Fig. 01 together with the UA1 and CDF 
measurements [5|, H| ■ The results are compared to 
a theoretical evaluation in leading-order pQCD [8j. 
The expectation is consistent with the data to within 
the combined experiment statistical and systematic 




Gluon Jet Energy (GeV) 

FIG. 4: Gluon splitting rate to charm pairs as a func- 
tion of the gluon jet energy. Measurements from STAR, 
UA1 [5( and CDF [6j collaborations are compared with 
pQCD calculations [3] using the indicated values of the 
charm quark mass (in GeV/c 2 ), Aqcd = 300 MeV, and 
a Peterson fragmentation function with e c — 0.06. 

uncertainties. Although the agreement is not strong, 
the conclusion that R g —, cS is small for energies ac- 
cessible at RHIC is clearly supported. The use of the 
PDG estimate for the fraction c — > D* + [l4| for the 
UA1 and CDF data does not change this conclusion. 

The pQCD expectation for i? g _ c5 [|[ can be com- 
bined with the STAR measured mid-rapidity jet dif- 
ferential cross section [l7|, covering jet pt in the 
range of 5 to 50 GeV/c, and the gluon jet frac- 
tion [22| to estimate the gluon splitting contribution 
to the charm production cross section. The gluon 
splitting contribution can in this way be determined 
for charm pt in the range of 2 to 10 GeV/c. The 
result is smaller by an order of magnitude than the 
pQCD expectation 23] for the total charm differ- 
ential production cross section at RHIC. It is thus 
small also compared to the measured total charm 
cross section [21(, which exceeds the pQCD expec- 
tation. 

In summary, we report the first measurement on 
the charm content in jets from p + p collisions at 
y/s = 200 GeV. The ratio N(D*+ + D*~)/N (jet) is 
measured to be 0.015 ± 0.008(stat) ± 0.007(sys) for 
D* mesons with fractional momenta 0.2 < z < 0.5 
in jets with a mean transverse energy of 11.5 GeV. 
This is consistent with perturbative QCD evaluation 
of gluon splitting into a pair of charm quarks and 
subsequent hadronization into D* mesons. The as- 
sociated cross section is smaller than the total charm 
production cross section at RHIC. We thus infer that 
the charm content in jets at RHIC energies has a 
small contribution from gluon splitting and is dom- 
inated by jets initiated by charm quarks. 
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